Sequential Strong Measurements and Heat Vision 
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We study scenarios where a finite set of non-demolition von-Neumann measurements are available. 
We note that, in some situations, repeated application of such measurements allows estimating 
an infinite number of parameters of the initial quantum state, and illustrate the point with a 
physical example. We then move on to study how the system under observation is perturbed after 
several rounds of projective measurements. While in the finite dimensional case the effect of this 
perturbation always saturates, there are some instances of infinite dimensional systems where such 
a perturbation is accumulative, and the act of retrieving information about the system increases 
its energy indefinitely (i.e., we have 'Heat Vision'). We analyze this effect and discuss a specific 
physical system with two dichotomic von-Neumann measurements where Heat Vision is expected 
to show. 

PACS numbers: 



Classically, the ability to extract information from a 
physical system is only constrained by our technological 
means: ideally, a Maxwell demon could have a record 
of all the relevant variables of a given experiment, that 
could be actualized regularly without disturbing the dy- 
namics of the system. In contrast, in quantum mechanics, 
measurement is an active process that inevitably modifies 
the system under observation. It is precisely the change 
the system faces during such a process that prevents us 
from measuring two non-commuting observables simul- 
taneously. 

Among the family of Positive Operator Valued Mea- 
sures (POVMs) used to model the measurement process 
in quantum mechanics, projective or von Neumann mea- 
surements stand out due to their special properties. Von 
Neumann measurements do not only retrieve informa- 
tion, they also establish a property on the system un- 
der observation: once a measurement is made and a re- 
sult obtained, any subsequent measurement will always 
output the previous value. Moreover, repeated measure- 
ments of the system will not modify its state any further. 
Von Neumann measurement are thus the closest analog 
of non-perturbative classical measurements, in the sense 
that they do not introduce external noise: any perturba- 
tion a system may experience is just a result of the state's 
non-definiteness of the property we want to estimate. 

This last feature, the fact that consecutive repetitions 
of von Neumann measurements do not alter the physi- 
cal state of the system, is the basis for the Zeno effect 
[l[ , where the evolution of a quantum system is literally 
frozen by a quasi-continuous measurement process. The 
Quantum Zeno effect has been shown recently to serve as 
a cooling mechanism in certain experimental situations 
■ On the other hand, from a tomographic point of view 
the Zeno effect is not very interesting, since the outcomes 
of all measurements performed are due to be the same. 

In this work, we are interested in the compromise be- 
tween the potential of von Neumann measurements for 
tomography and the perturbation of the quantum sys- 
tem over which they are implemented. We have just 
seen that allowing only one type of measurement leads 



to a very poor dynamics and not-very-informative to- 
mography. This motivates us to study physical scenarios 
where a limited (finite) number of projective measure- 
ments with a finite number of outcomes each are avail- 
able. 

Let us first explore how the fact that we have several 
measurements to choose from can suppose a huge differ- 
ence when we want to do tomography. 

Quantum tomography 

In a classical system, the ability to perform two differ- 
ent measurements of d outcomes each can only retrieve, 
at most, 21og 2 (d) bits of information. Indeed, if we mea- 
sure the height of an individual and then its weight, we 
shall not expect to gather more information by measur- 
ing its height one more time. However, in the quantum 
world, two measurements +, — will not commute in gen- 
eral. Given three instants of time to < t\ < t2, this can 
lead to 'paradoxical' situations where the output of mea- 
surement + at time to may differ from the output at time 
ti if there has been an intermediate measurement of — 
at time t\. 

Sequential measurement schemes have already proven 
useful in the weak measurement scenario [H, and could 
play an important role in future experimental tests of 
Quantum Mechanics Q. However, the possible appli- 
cations of consecutive 'strong' non-commuting measure- 
ments for tomography seem to be absent from scien- 
tific literature. This is certainly peculiar if one takes 
into account that, giving two non-commuting projectors 
P,P', the operator space spanned by the POVM el- 
ements {P,PP'P,P'PP',PP'PP'P,...} describing con- 
secutive measurements may well be infinite dimensional, 
and so the corresponding estimated probabilities can give 
us access to an infinite number of parameters character- 
izing the state of our system. The current experimen- 
tal credo, though, does not echo this simple observation: 
even when non-demolition interactions are available, sys- 
tems are typically discarded after a single-shot measure- 
ment (e.g., in Wigner function estimation 5]). Exploring 
many non-commuting properties of an unknown quan- 
tum state thus involves an unnecessarily large number of 
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system initializations. 

The previous ideas can be better illustrated by refer- 
ring to a specific physical system. 

Consider a spin 1/2 neutral particle in a rectangular 
box of sizes L x Lq x Lq, with one of its vertices situated 
at the origin of coordinates. We will assume that Lq <C 
L, so this system can be regarded as a one-dimensional 
object subject to the potential 



V{x) 



0, for < x < L, 
oo, otherwise. 



(1) 



The hamiltonian describing the evolution of such a par- 
ticle is thus 



H = 



2m 



V(x) 



(2) 



where m denotes the mass of the particle. The eigenvec- 
tors of this system are {|n) : n € N + } with (x\n) =: 

^n(x) — y^sin( ZL j £ ), each with associated energy 

E n = 7r 2 n 2 /2mL 2 . The spin degree of freedom of the 
particle can be modeled through a 2-dimcnsional Hilbert 
space C 2 . Along these pages, the three pairs of kets 
{|0),|1)}, {|+),|-)} and will represent the 

eigenstates of the Pauli matrices o~ x , o~ y and er z , respec- 
tively. It follows that the Hilbert space T-L <S> C 2 where 
this particle lives can be then expanded in the basis 
{|n, ±i) : n e N + }. As usual, quantum states in this sce- 
nario will be regarded as positive semidefinite normalized 
elements of S\ = {A : tr(|A|) < oo}, the set of trace-class 
operators. Not to be confused with the class of opera- 
tors S2 = {A : tT(AA*) < 00}, that will also play an 
important role shortly. 

Suppose now that our technology allows performing 
almost instantaneous von Neumann spin measurements 
along the u z direction over such a particle without affect- 
ing its canonical degrees of freedom. Defining \<t> ± (x)) = 
cos(fca;)|0)±sin(fcx)|l), then, the two von Neumann mea- 
surements + and — associated to the projectors 



F ± = / |a:)(a:|® ^(aOX^aOlda: 



(3) 



can be physically realized by applying a magnetic field 
±B along the u y direction with an intensity varying 
linearly with the x coordinate }14| . measuring the spin 
and then applying the inverse field =p_B. Indeed, if 
B = —(bx)uy and the magnetic interaction is me- 

y,z=0 

diated through a hamiltonian H s 
check that 



-pB ■ <7, one can 



a ^iH a At 



(I®|0}<0|) 



+iH s At 



(4) 



provided that we switch the field B for a time At = -K. 
This time will have to be very short (and so the magnetic 



density b will have to be very strong) if we want to neglect 
the evolution of the system due to the main term @ 
during At. 

Module technicalities, we are thus able to implement 
two different von Neumann measurements over our sys- 
tem. We will call the outputs of such measurements 
and 1 when the related projectors are F ± and I - F ± , 
respectively. 

We will next show that the statistical analysis of re- 
peated measurements allows reconstructing the proba- 
bility density p(x) of the particle inside the box. Before 
proceeding, note that we are not making any assumption 
on the initial coupling between the spin of the particle 
and its canonical degree of freedom (i.e., they could be 
classically correlated, or even entangled). 

Suppose that we measure + and — alternatively N 
times. For any a € {0, 1} N , define the probabily P(a\ + 
— I — ...) of obtaining outcomes (01,02,03,...) during a 
sequential performance of the measurements +,—,+,... 
Then it is easy to check that 



P(0,0,0, 



+ ...) + P(l,l,l,. 

<cos 2 ^ 1 )(2fcx)). 



(5) 



Using the fact that cos 2M (#) is a linear combination of 
the functions {cos(2j6>) : j = 0, ...,M}, we can therefore 
estimate the values {(cos(4fcjx)) : j = 0, N — 1} from 
the statistical analysis of N repeated measurements. If 
k < we can then extend p(x) to an even function de- 
fined in [— 7r/4fc, ir/ik] and use our statistical knowledge 
to infer the coefficients {cj}j of the Fourier expansion 



p(x) 



Cj cos(4fcjx), for x € [0, L]. 



(0) 



At first glance, one could argue that, since the proba- 
bilities appearing in eq. (J5J) decrease exponentially faster 
with N, it would require a large number of samples to es- 
timate them. However, a proper calculation shows that, 
as long as p(x) has a finite slope at x = 0, the right-hand 
side of eq. (0 decreases as 0(1/N). 

We have just proven that, as we repeat measurements 
+ and — , we obtain more and more useful information 
about our initial state. However, as we will see soon, this 
information comes at a price. 

Heat Vision 

As we said at the beginning, we are interested in how 
the state of the system is affected by the measurement 
process. Of course, one would expect such a state to 
depend on the implemented measurement strategy, the 
process by which we choose which measurement to apply 
at a given time. 

Any scenario where a limited number of von Neumann 
measurements are available can be studied by analyz- 
ing the statistics that result from randomly applying 
one measurement or the other. It is thus legitimate to 
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consider the behavior of the system under independent 
random measurement strategies, where the probability 
p x > of performing measurement x is the same on 
each round. We will see that the overall effect of random 
measurement strategies can have a very different nature 
depending on whether the dimension of the underlying 
Hilbert space is finite or infinite. 

Suppose, indeed, that we perform with probability p x 
a measurement x G {1,2, ...,s} defined by the complete 
set of projectors {F£ : a — 1,2, ...,d} C B(H), for some 
Hilbert space H. The action of the resulting map d over 
an initial state a would then be given by 

Sl(a) = Y J P x F>K- (7) 

a;, a 

Note that we can see a = 53ijK)0'l ^ ^2 as an element of 
T-L®T-L via the isomorphism [15j a — >• \a) = JV^. &ij K) \j) ■ 
d can then be regarded as a superoperator d £ B('H'^'H), 
given by 

n = (8) 

x,a 

Seen as an operator, d is both hermitian and positive 
semidefinite. Moreover, since this map is also unital, 
its operator norm will be upperbounded by 1 [fj|. It 
follows that the spectrum of d is in [0, 1] and so the 

— N 

limit lim^v-i-oo " exists and is equal to IIi, the pro- 
jector onto the space of eigenvectors of d with eigen- 
value 1. Note that such a limit does not depend on 
the initial probabilities p x , as long as all of them are 
strictly positive. This is so because IIi is just the (pos- 
sibly null) projector onto the intersection of the spaces 

H x = spaS{|0) : Ea F a® m*\4>) = \m- 

Since in finite dimensions S2 and the set of trace-class 
operators coincide, we can conclude that repeated appli- 
cations of the mapping d will bring any quantum state 
a to a limiting state Sl°°(a) E Si. That is, even though 
the system may experience some perturbations at the be- 
ginning of the measurement process, given some time it 
will stabilize into a steady state. The finite dimensional 
case is thus resemblant of the quantum Zeno effect, i.e., 
repeated applications of an independent random mea- 
surement strategy will 'freeze' the state with respect to 
some types of environmental noise. 

In infinite dimensional systems, though, the norms || • || 1 
and || • || 2 are not equivalent, so the convergence of the 

sequence of vectors (d N \a)) in % <£> Ti does not neces- 
sarily imply the convergence of the sequence of states 

(d N (a)) in Si. In particular, it could happen that 
JV 

liniAr-^oc d \a) = G 82- 

As shown in Appendix |A1 if (d N (a)) does not con- 
verge in Si, such a sequence of quantum states will not 
admit a finite dimensional approximation. Moreover, the 
sequence of mean values tr(0 JV (a)H) can be shown to di- 
verge for any 0-band hamiltonian H iff d N (a) does not 



converge in Si. In such a situation, the action of gath- 
ering information about the system will then increase its 
energy up to arbitrarily high values. We will hence call 
such an effect 'Heat Vision', in analogy with Superman's 
famous ability to induce heat with the power of his stare 
0- 

We will next prove that, remarkably, our previous ex- 
ample exhibits Heat Vision for any initial state a G Si 
we place as an input. 

First, since Heat Vision does not depend on the actual 
weights we assign to each measurement, we can assume 
w.l.o.g. that we perform any of the two measurements 
with probability 1/2. The measurement channel d is thus 
equal to 

(1(a) = i{ J2 F3(jF3 + ( f ~ Fi ) a $ ~ F3 ^- ( 9 ) 

j=+- 

As shown in Appendix [B] for any state a, 

lim tr{[fi*(<r)] a } = 0. ( 10 ) 

N— >oo 

Clearly, the values (\\(l N (a) — 0||i) do not tend to 0, and 
so the sequence of states (d N (a)) does not converge in 
trace norm. Also, equation (TT0|) . together with the Renyi 
inequality S(x) > — log2(tr{x 2 }), valid for any normal- 
ized quantum state, implies that the von Neumann en- 
tropy of (l N (a) tends to infinity. 

As we mentioned earlier, the fact that (d N (a)) does 
not converge in Si implies that the energy of the system 
will diverge for any input state, as long as the hamilto- 
nian describing the system has no energy bands in its 
spectrum. Such is the case of the hamiltonian described 
by eq. @, and, indeed, one can check that the kinetic 
energy of a quantum state after N applications of the 
channel d obeys the formula 

EW = + —N. (11) 

m 

That is, on average, the temperature of the system will 
grow linearly with the number of measurements per- 
formed. 

Let us briefly recapitulate what is happening here: in 
order to perform measurement ±, first we have to apply 
a strong magnetic field ±B for a time. It is not surpris- 
ing then that the energy of the system increases when we 
perform such a change. However, after the spin measure- 
ment we will apply the opposite field, ^fB, thus inverting 
the previous unitary process. The only reason why the 
energy of the system (or for the same sake, the state of 
the system) changes is thus that we are measuring a sin- 
gle qubit in between! Indeed, the system is engineered 
in such a way that the entropy associated to such a spin 
measurement accumulates and accumulates in the canon- 
ical degree of freedom of our particle until the setup can- 
not stand more energy. This behavior is to be compared 
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FIG. 1: Dynamics of a quantum system subject to repeated 
von Neumann measurements. 

to the finite dimensional case, where it is always possible 
to find states of entropy 1 that are invariant under the 
action of any two dichotomic measurements jL6 ( ]. 

But Heat Vision can manifest in even more extreme 
ways: in Appendix [D] we describe systems with 5 di- 
chotomic observables where the purity of any initial state 
a subject to N sequential measurements is bounded by 
X 2N , where A < 1 is independent of a. The extropy of 
such states will thus increase at least linearly in N. The 
measurements involved, though, are quite abstract, and 
most likely impossible to implement in any present lab- 
oratory. 

By showing that Heat Vision does emerge (at least, 
theoretically) in some physical systems, now we have a 
clear picture of how quantum systems evolve through se- 
quential von Neumann measurements. The conclusions 
are illustrated in Figure [TJ where the finite and infinite 
dimensional case are differentiated. 

It is tempting to think that Heat Vision arises in some 
infinite dimensional systems and not in finite dimensional 
ones just because in the former we can extract an infinite 
amount of information about the initial state of the sys- 
tem. However, in Appendix [Cl we describe a system with 
two dichotomic von Neumann measurements where the 
statistical analysis of sequential observations retrieves in- 
finitely many parameters, but Heat Vision never shows 
up. Heat Vision is thus not equivalent to the possibility 



of accessing an infinite amount of information, but an 
independent property of the system under study. 

Finally, we would like to point out that both tomogra- 
phy with strong sequential measurements and Heat Vi- 
sion can be experienced with current technology. Indeed, 
consider a regular ion trap setup with just one ion. Then 
two internal states of such an ion could play the role of 
the spin in the system described above, while the dis- 
placement of the ion along the trap could account for 
the canonical degree of freedom. 'Spin' measurements 
in this scenario can be performed in the standard way, 
i.e., exciting one of the internal levels with a laser and 
counting the number of emitted photons. Analogously, 
an interaction of the form H s can be induced by a laser 
beam in standing wave configuration [9]. In a usual ion 
trap setup, the former dispositions would effectively im- 
plement the measurements Q over a particle subject to 
a harmonic potential; in order to recreate the square po- 
tential, a trap of the form [l(| can be used. Due to the 
Brans-Dicke approximation, though, the whole experi- 
ment must be conducted in the low temperature regime 
(up to - IK). 

Conclusion 

In this paper we have studied the use and effect of re- 
peated von Neumann measurements. We have pointed 
out an extreme scenario where the statistical analysis of 
binary outcomes of sequential measurements allows es- 
timating the full probability density of a trapped parti- 
cle. We have also shown how sometimes the action of 
alternating measurements can lead to a non-convergent 
dynamics in infinite quantum dimensional systems. This 
phenomenon always comes together with an unbounded 
energy increase (hence its name 'Heat Vision'), and can 
be observed experimentally in current ion trap setups. Is 
this the end of the story? Dreaming on, one could con- 
ceive a new architecure for quantum computing based 
on sequential strong measurements. In this model, the 
user could perform a (small) number of non-commuting 
dichotomic measurements over a continuous variable sys- 
tem, and computations would be carried out by deciding 
which measurement to perform at every step. Although 
still vague, we hope to explore this idea in future com- 
munications. 

The authors would like to thank A. Retzker for sug- 
gesting experimental implementations. This work was 
supported by the Spanish grants I-MATH, MTM2008- 
01366, S2009/ESP-1594 and the European projects 
QUEVADIS. 



[1] E. C. G. Sudarshan, and B. Misra, Journal of Mathemat- 
ical Physics 18, 4, 756 (1977). 

[2] N. Erez, G. Gordon, M. Nest, G. Kurizki , Nature 452, 
724-727 (10 April 2008) 

[3] G. Mitchison, R. Jozsa, and S. Popescu, Phys. Rev. A 



76, 062105 (2007). 

[4] T. Fritz, J. Math. Phys. 51, 052103 (2010). 

[5] L. G. Lutterbach and L. Davidovich, Phys. Rev. Lett. 78, 
2547 (1997); G. Nogues, A. Rauschenbeutel, S. Osnaghi, 
P. Bertet, M. Brune, J. M. Raimond, S. Haroche, L. G. 



Lutterbach and L. Davidovich, Phys. Rev. A 62, 054101 
(2000). 

[6] D. Perez-Garcia, M. M. Wolf, D. Petz, M. B. Ruskai, J. 

Math. Phys. 47, 083506 (2006). 
[7] See, for instance, Burn!, Superman, 218, DC comics 

(2005). 

[8] J. Clerk-Maxwell, Treatise on Electricity and Magnetism, 
page 319, Oxford: The Clarendon Press. ISBN 0-486- 
60636-8 (1873). 
[9] D. F. V. James, Appl. Phys. B 66, 181 (1998). 
[10] P. Staanum, S. D. Kraft, J. Lange, R. Wester, and M. 

Weidemller, Phys. Rev. Lett. 96, 023201 (2006). 
[11] R. Renner. Security of Quantum Key Distribution, PhD 
thesis, Swiss Federal Institute of Technology (ETH) 
Zurich (2005). 

[12] D.V. Voiculescu, K.J. Dykema, A. Nica, Free random 
variables. A noncommutative probability approach to free 
products with applications to random matrices, opera- 
tor algebras and harmonic analysis on free groups, CRM 
Monograph Series, 1, American Mathematical Society 
(1992). 

[13] G. Pisier, An Introduction to Operator Spaces, London 
Math. Soc. Lecture Notes Series 294, Cambridge Univer- 
sity Press, Cambridge (2003). 

[14] This can be implemented by creating a constant gradi- 
ent field of the form B — —byu x — bxu y by means of a 
Maxwell coil Q. Since the particle's Y coordinate is null, 
the effective magnetic field will be B = —bxu y . 

[15] Indeed, notice that tr(cr'o") = (o~\a). 

[16] In a finite dimensional system, any pair of projectors F + 
and F~ can be simultaneously 2 x 2-block-diagonalized, 
i.e., there exists a basis where F + = ® n F„ , F~ — ©„F~, 
with F+,F~ € M 2X 2 (Jordan, 1875). Any state of the 
forma = 2X 2ffi...©y ©•••©0 2X 2 thus satisfies F ± oF ± + 
(I-F ± )o-(I-F ± ) = g. 



6 



Appendix A: Convergence in Si 

Theorem 1. Let <j( N ^ be a sequence of normalized quan- 
tum states. Then, the following conditions are equivalent: 

1. limjv^oo exists in Si. 

2. For some (and thus for all) arbitrary countable or- 
thonormal basis {|^)} ofH, the limits 

c n , m = lim tr{o-W\n)(m\} (Al) 
exist and are such that 



n=0 

Proof. Suppose that 1) is true. Then, there exits an el- 
ement a £ S± such that limjv^oo \\& — cr( N '\\i = 0. We 
remind the reader that, for any self-adjoint operator A, 

\\A\\i= sup tv{A-X}. (A3) 

I>X>-I 

Now, let {|n}} be any orthonormal basis for H. The 
operator \n)(n\ satisfies I > ±|n)(n| > —I, so 

II* - °- (N) \\i > \tr{\n)(n\a} - tr{\n)(n\a^}\. (A4) 

It follows that the limit limjy->-oo tr{erW|n)(n|} exists 
and is equal to tr{<r|n)(n|}. Analogously, from the re- 
lations 

I > |ra)(m| + \m)(n\ > -I, 

I> i(\n){m\ - \m){n\) > -I, (A5) 

it can be shown that limjv->oo tr{<rW|n)(m|} exists as 
well. Finally, I > ±1 > —I, which means that \\<r — 
<J {N) \\i > \tr{a-trW}\. Since tr(oW) = 1,VJV, we 
have that tr(o-) = 1, and so J2 n c n , n = tr{^ n \n)(n\a} = 
tr(cr) = 1, and 2) is proven true. 

Conversely, suppose that 2) is true, and consider the 
operator 



o = ^ c n,m\m)(n\. (A6) 

n,m— 

This operator is bounded. Indeed, let G 
span{|n)}. Then, 

\{v\a\w)\ = lim \ti{o- (N) \w)(v\}\ < J(w\w)(v\v). 

(A7) 



Likewise, it can be shown that (v|cr|v) > for all 
\v) G span{|n)}, i.e., a > 0. Moreover, by equation (IA2|) 
tr{a} = 1, so a € Si. 

Let Pk = 2 n= o l n X n l- Then equation (|A2p im- 
plies that, for any e > 0, there exist K,M such that 
tr( (7 ( JV )p if ) > 1 - e.ViV > M. Applying twice the rela- 
tion [HI 

\\p-PkpPk\\i < 2Vtr( i o){tr( ( o) - tr(pP^)}, (A8) 
valid for any p > € Si , we have that 

lk (JV) -*||i <^+\\PK^P K -P K aP K \\ u (A9) 

for N > M. Note that the last term of Eq. (|X9j) tends to 
as TV tends to infinity (because we are evaluating the 
trace distance between two K + 1 x K + 1 matrices that 
converge entry- wise). It follows that lim n ->oo H"' — 
^ 4ye. Since e was arbitrary, we conclude that 
lini/v->oo ll '^ — *||i = 0, and so (cr^) converges in 

□ 

In the particular case where erW = Q N (a), for some 

initial state a and some channel f2 with I > Q > 0, 

TV 

the existence of the limits (jAll) is automatic, since 
converges in B(H <g> "H) and |n)|m)* £ H ® H. This 
implies that convergence in Si in that case is equiva- 
lent to the existence of a basis {|n) : n £ N} such that 
limx^oo limjv-i-oo tr(Pif 0^(0-)) = 1 (note the order of 
the limits). If the latter is the case, then can always 
be described by a finite dimensional system, i.e., for suf- 
ficiently high K, we can approximate by the state 
P k o- (N) Pk G B(C k+1 ), for all N. 

In order to establish a connection between energy and 
convergence in Si in realistic scenarios, we will have to 
restrict the usual definition of hamiltonian. 

Definition 2. O-band energy operator 

Let E be a self-adjoint operator acting over an infinite 
dimensional (separable) Hilbert space %. We will say 
that E is a 0-band energy operator iff 

1. The spectrum of E is discrete. 

2. For any E € K, there is only a finite number of 
linearly independent eigenvectors of E with eigen- 
values less or equal than E. 

Examples of 0-band energy operators are the harmonic 
oscillator, the double- well potential, a particle in a box... 
and, more generally, the hamiltonian of any finite number 
of particles subject to a potential that can be bounded 
from below by a harmonic trap. Since most experiments 
in Physics are not performed in the open air, but inside 
closed chambers, 0-band energy operators provide a very 
good dynamical description of those quantum systems 
accessible in the lab. 

The next lemma relates convergence in Si with energy 
considerations. 
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Lemma 3. Let (crW) be a sequence of normalized quan- 
tum states such that limjv^oo (T exists in S^. Then, 
linijv-^cx) crW exists in S\ iff, for some O-band energy 
operator E Yrni^^oa tr{Ea^ N ^} exists. 

Proof. Suppose that the sequence (cr^) does not con- 
verge in Si , and let E be an arbitrary O-band energy oper- 
ator E — X^^Lo E n \n)(n\, where {|n) : n e N} are a basis 
of eigenvectors of E, with energies Eq < Ex < E2 < .... 
Let Pk be the projector Pk = J2n=o l n X n l- We will 
prove that there exists 1 > A > such that, for any 
E > E , liniAT^oo tr{£JcrW} > (1 - X)_E a + XE. 

Indeed, let K be such that Ek+i > E. If liniAr^oo ffW 
does not converge in Si , by Theorem [T] we have that 
limjv^oc tr{a^ N) P K } < lim^oo liniAr^oo tr{cr (A,) P„} 
1 A. with A > independent of K. On the other hand, 



If lim^ K(s) = K < 00, then E„=oP W (™) = 1 
for all N. We can thus simply define the O-band energy 
operator E = Xri=if+i n l n )( n l an( ^ we wou ld have that 
tv{a^E} = < 00. 

Suppose, on the contrary, that lim^oo K(s) = 00, and 
define the sets of natural numbers Jo = [0, 1^(1)] and 

I s = [K(s) + l,K(s + l)l if K (s) + 1 < K(s + 1) 

0, otherwise, (A14) 

for s > 1. These sets are finite and disjoint, and satisfy 
U^Lo^s = N. Denoting the projector EneJ \ n )( n \ as 
we thus have that the positive operator 



£ = J2 ^ Sp * ( A15 ) 

s=0 



£ = ^£ n |n)(n| >£ 0j Pk + £(I-Pk), (A10) 
so we arrive at 



is 0-banded. 
Finally, notice that 



lim tr{£cr (Ar) } > (1 - X)E + XE. (All) 

JV— »oo 

Conversely, suppose that (o-^ N ') converges to the nor- 
malized state a. We will prove that there exists a 0- 
band energy operator E such that tr{a^ N 'E} is bounded 
and linijv->oo tr{a^ N 'E} exists. Let {\n) : n e N} 
be any basis for % and define the probability distribu- 
tions p( N \n) = (n|<rW|n), p(n) = (n\a\n). Also, let 
K : N —> N be the mapping defined as 



K 1 
K{s) = {minif > : ^p W (n) > 1-— ,ViV}. (A12) 

It is important to note that K(s) < 00 for all s. Indeed, 
suppose that K(s) = 00, for some s. That would imply 
that, for any number K, there exists an TV such that 
En=oP (J ^( n ) <1 ~W- Now > for < e < 1/2 S+1 , choose 
M such that \\a - cr (Ar) || 1 < e, for all N > M, and choose 
K such that 



K 1 
^ p W(„)>l-_ + e> VJV<M 1 

n=0 



K 

E 

71 = 



P» > 1 - ^ + 2e. 



(A13) 



(JV) 



Then, for N > M, | E*= P (JV) (») - P»l < Ik 

< e. It follows that the expression Eri=o ( n ) > 
1 — ^ + e holds for all AT, thus contradicting our initial 
claim. 

We will differentiate two cases depending on the exis- 
tence of the limit lim^oo K(s). 



EP W W< E P (N) (n)<l-, (A16) 

n=K(s) + l 



nel 



for s > 1. This implies that, for any A^, 



#)Etr(a^)<l + ^|j =1+ Vf L T' 

(A17) 

i.e., the sequence (E^ N >) is bounded. 

Now, define Q t = E s =o-^s- From all the above, it is 
clear that 



0<EW-tr(Q t *WQ t E)< £ . (A18) 

Likewise, 



< E - tr(Q t &Q t E) < f^j 

s=t+l V / 



where J5 = tr(o-.E'). It follows that 



(A19) 



\E - EW\ < \tr(Q t [& - o-W]Q t E)\ + (A20) 

Taking the limit N — > 00 and then t — ¥ 00, the 
right-hand-side of equation (|A20|) vanishes, and thus 
limiv->ootr(<r<W.E) = tv(aE). 

□ 

Let us make a final remark. 
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Lemma 4. Suppose that the channel I > f2 > 
has the property that, for any quantum state 
a, the sequence (fi (o~)) does not converge in 
S\. Then, limN-^oo fi (<t) = in S2, i.e., 
liixLTv-s-oo tr{(j( N ) |</))('(/'|} = for any pair of states 

\<t>)M)- 

Proof. Suppose that there exists a basis {\n) : n e N} 
for the Hilbert space such that the coefficients c„. m = 
1™^^^ tr (fl N (a)\n)(m\) do not satisfy condition (|A2|) . 
Then, following the proof of Theorem [TJ one could build 
an operator a > € Si such that tr(<r) < 1. Now, 
if tr(cr) 7^ 0, then a' = a/tr(a) would be a quantum 
state such that Q,(a') = a', contradicting the main as- 
sumption. We thus have that tr(er) = which, together 
with a > implies that a — (and so, (4>\a\tp) = 
lim A r^ oo tr{aW|V)(0|} = O, for all |V>>, □ 

Appendix B: Derivation of Equation (|10|) 

Viewed as a superoperator, the channel © can be seen 
equal to 

Tt = ^(1®U)J dxdy\x)(x\ ® \y)(y\ <g> 

<g>(M(x-y)® M(x + y))(JL®rf). (Bl) 

where 

„ r , . ( 1 cos[2fca;l \ 
M & = { cos[2 fca; ] 1 J 

= cos 2 [fc.T]|+X+l + sin 2 [fear] I -)(- 1, (B2) 

and 

U= |+t,-*){Q,0| + |-t,+*){0,l| + 

+|+* 1 +t)(l,0| + |-*,-t)(l,l|. (B3) 

It thus follows that 

— N 1 f 

n = -(1(g) U) I dxdy\x){x\ <g> \y)(y\ <g> 

<g> (M(x - y) N ® M(x + y) N ) (I <g> E/ t ),(B4) 
and one can then check that 

tr{[r! A '(K±i)(n,± l |)] 2 } = 

= (n, ±i, n, =pi|r2 2JV |7T,, ±i, n, ^fi) = 

= \ [ [ dxdy^ n {x) 2 ^ n {yf ■ 
z Jo Jo 

• (cos 4JV [k(x - y)} + sm 4N [k(x - y)]) . (B5) 



Taking into account that n {x) 2 ^! „(y) 2 < 4/L 2 and 
changing to variables x' = (x — y)/L,y' = (x + y)/L, 
we have that 

tr{[n N (\n,±i){n,±i\)] 2 } < ip{N) 2 := 
= 2 - \x\] ■ [cos iN {kLx)+sm 4N (kLx)]dx.(B6) 

For any initial state cr, applying the Schwartz inequality 
one arrives at 

tr{n N (a)\n,±i)(n,±i\} = tv{aft N (\n, ±i)(n, ±i\)} < 

< y/i^yJtT{[n N (\n,±i)(n,±i\)] 2 } < <p(N), (B7) 

with limjv-»oo f(N) = 0. Note that this expression is 
independent of n: the occupation of each of the states 
{|n,±i)} tends uniformly to zero. In other words, the 
energy density distribution of a quantum state subject 
to sequential measurements + and — neither converges 
nor displaces, but flattens. 

Now, Theorem [T] in Appendix [5] states that any se- 
quence (Q N (a)) does not converge in Si iff there ex- 
ists some orthonormal basis {\ip n ) ■ n & N} such that 
J2 n ^N^oo(^n\^ N (<j)\ipn) < 1. Moreover, by lemma 0] 
in the same Appendix, if such is the case for any initial 
state a, then limjv->oo £l N (o~) tends to in the || • H2 norm. 
It thus follows that 

lim tr{p}*(<r)] 2 } = 0, (B8) 

N— >oo 

for all initial states a. 



Appendix C: An example of unbounded extraction 
of information with no heat vision 

Consider a separable Hilbert space H, and let {\n) : 
n G N} be an orthonormal basis for %. Then, we can 
define projector operators acting over W®C 2 as 

00 

G ± = J2\n)(n\^\^ ± (n))( ( f > ± (n)\. (CI) 

71=1 

The analog of equation §5§ follows straightforwardly, and 
so we can combine our projective measurements to esti- 
mate the mean values {(cos(4fcjn)) : j £ N}. Choos- 
ing k irrational, the statistical analysis of repeated mea- 
surements can thus provide us with an infinite amount 
of information about the ocupation number distribution. 
However, let fl be the channel that results when we ran- 
domly apply one measurement or the other. Then, for 
any initial quantum state <r, it can be shown that 

00 

lim n N (a) = y^(n\tt C 2(a)\n)\n){n\ ®I 2 /2, (C2) 

N— too ^ — ' 

n=l 
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that is, the system does not exhibit heat vision in any 
case. Moreover, if the energy operator is diagonal in the 
{\n) : n £ N} basis, the energy of the system does not 
even vary during the measurement process. 

Appendix D: Extreme Cases of Heat Vision 

An extreme case of heat vision can be found in the 
following system. Consider the group G that results out 
of the free product [l2[ of Z2 with itself s times, i.e., 
s times 

G = Z2 * Z2 * ... * Z2, and take /2(G) to be our Hilbert 
space. 

In this context, the left regular representation of an 
element g £ G is defined as the unitary operator X(g) : 

hip) -> h(G) such that A (5)l.9'> = Iffff')) for an y 9' e G 

jl3{ . The left regular representation of the generators 
A(<?i) thus satisfies X(gi) 2 = 1, A(<fr) = \(gi)^. This im- 
plies that, for each generator gi, there exists an asso- 
ciated projector (A(<fo) + I)/2, and so each gi defines a 
quantum dichotomic measurement. The channel VL that 
results when we perform one of the s measurements with 
probability l/s can be written as 

11 =\ (jE A ^)® A ^)* +I ) • ( D1 ) 

Define the operator IT as the projection onto the sub- 
space of h(G) spanned by all the elements of G that start 
with the symbol gi, and note that X(gi) = x% + yi, with 
Xi = X(g l )\l l ,y l = UiX(gi). We have that 

|lX>Gfc)®A(ff 4 )l = ||^A(^)|| < ||£>|| + ||X>||, 

i—1 i—1 i—l i—1 

(D2) 



where in the first equality we have made use of Fell's 
absortion principle [13]. On the other hand, for any two 
sets of operators {Ai} {Bi}, 

iiE^ii ^ iiE a ^ii 1/2 iiE b 1^h 1/2 - (° 3 ) 

i i i 

Taking (A; = I, B t = Xi) and (Ai — yi, Bi = I), we have 
that the last term of equation (ID2[) is upperbounded by 



Vi(jiE n # /2 + iiE n *ii 1/2 J' (° 4 ) 

that, in turn, is upperbounded by 2y/s, since J^. IT < I. 
It follows that 



||fi|| < 1/2 + 1/V5. (D5) 
The norm of fl as an operator in h(G) ® /2(G) is thus 
smaller than 1 whenever the number of measurements 
is greater than 4. This phenomenon can only occur 
in infinite dimensional systems, since, for any finite di- 
mensional unital map uJ on C d <S> C d , S7|I C ;) = |IA, for 

\ J d) = Ei=l and SO = !• 

Suppose, then, that s > 5 and so ||0|| = A < 1, and 
let a £ Si (£2 (GO) be any arbitrary normalized quantum 
state, with <j( N > = Q N (a). Our previous discussion im- 
plies that 

tr{{a^} 2 } < X 2N tv{[a} 2 } < X 2N . (D6) 

That is, the purity of any initial state decreases exponen- 
tially with the number of applications of the channel, and 
so the system exhibits Heat Vision for any input state. 



